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Abstract 

We introduce the new concept of "metaiines" for manipulating the amplitude and phase 
profile of an incident wave locally and independently. Thanks to the highly confined graphene 
plasmons, a transmit-array of graphene-based metaiines is used to realize analog computing 
on an ultra-compact, integrable and planar platform. By employing the general concepts of 
spatial Fourier transformation, a well-designed structure of such meta-transmit-array combined 
with graded index lenses can perform two mathematical operations; i.e. differentiation and 
integration, with high efficiency. The presented configuration is about 60 times shorter than 
the recent structure proposed by Silva et al. (Science, 2014, 343, 160-163); moreover, our 
simulated output responses are in more agreement with the desired analytic results. These 
findings may lead to remarkable achievements in light-based plasmonic signal processors at 
nanoscale instead of their bulky conventional dielectric lens-based counterparts. 

Recently, realization of analog computing has been achieved by manipulating continuous val¬ 
ues of phase and amplitude of the transmitted and reflected waves by means of artificial engineered 
materials, known as metamaterials, and planar easy-to-fabricate metamaterials with periodic arrays 
of scaterrers, known as metasurf aces.Both above-mentioned platforms offer the possibility of 
miniaturized wave-based computing systems that are several orders of magnitude thinner than con¬ 
ventional bulky lens-based optical processors.^ 
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Challenges associated with the complex fabrication of metamaterialsl^ besides absorption 
loss of the metal constituent of metasurfaces,degrade the quality of practical applications of 
relevant devices. As a result, graphene plasmonics can be a promising alternative due to the tunable 
conductivity of graphene and highly confined surface waves on graphene, the so-called graphene 
plasmons (GPs).^^^^^ 

We present a planar graphene-based configuration for manipulating GP waves to perform de¬ 
sired mathematical operations at nanoscale (see Figure [T]). By applying appropriate external gate 
voltage and a well-designed ground plane thickness profile beneath the dielectric spacer holding the 
graphene layer, desired surface conductivity values are achieved at different segments of graphene 
layer. 115122 ] illustrate the applications of the proposed configuration, two analog operators, i.e. 
differentiator and integrator, are designed and realized. The proposed structure will be two dimen¬ 
sional which is an advantage compared to the previously reported three dimensional structures'^ 
that manipulate one dimensional variable functions. 

We introduce a new class of meta-transmit-arrays (MTA) on graphene: metalines which are 
one dimensional counterpart of metasurfaces. Our approach for realizing mathematical operators 
is similar to the first approach of [^; i.e. metaline building blocks (instead of metasurfaces) per¬ 
form mathematical operations in the spatial Fourier domain. However, the main advantage of our 
structure is that it is ultra-compact (its length is about 1/60 of the free space wavelength, X) in 
comparison with the structure proposed by Silva et al. whose length is about A/3.^It should be 
also noted that, in this work, the whole structure (including lenses and metalines) is implemented 
on graphene. Therefore, the total length of the proposed device is about A/4, about 60 times 
shorter than the device reported in [^. 

THEORETICAL ERAMEWORK 

The general concept of performing mathematical operation in the spatial Fourier domain is graph¬ 
ically shown in Figure]^ In this figure, z is the propagation direction, h{x,y) indicates the desired 
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two dimensional impulse response, f{x,y) is an arbitrary input function and g{x,y) describes the 
corresponding output function. The whole system is assumed to be linear transversely invariant 
and thus the input and output functions are related to each other via the linear convolution:^^ 

g{x,y) = h{x,y)*f{x,y) = j J h{x-x',y-y')f{x',y')dx'dy' (1) 

Transforming eq ?? to the spatial Fourier domain leads to: 

G{k,,ky)=H{k,,ky)F{k,,ky) (2) 

where G{kx,ky), H{kx,ky) and F{kx,ky) are the Fourier transform of their counterparts in eq ??, 
respectively, and {kx,ky) denotes the 2D spatial Fourier domain variables. 

For our two dimensional graphene-based system, as shown in Figure]^ the above formulation 
is sufficient to be presented in one dimension. In this figure f{x) and g{x) represent transverse field 
distribution of the incident and transmitted waves, and H{kx) is the appropriate transfer function. 
Accordingly, eq ?? can be interpreted as: 

g(x)=IFT{//(fc,)FT[/(v)]} (3) 

where (I)FT means (inverse) Fourier transform. The transfer function is given in the spatial 
Fourier domain kx, and, on the other hand, the incident wave {f{x)) is also transformed into 
the Fourier domain. Hence, any transfer function can be realized by properly manipulating the 
Fourier-transformed wave in its transverse direction 

Fourier transform can be carried out by a lens at its focal point. On the other hand, since 
realizing inverse Fourier transform with real materials is not possible, so instead of eq ??, the 
following relation should be used: 


g(-x)=FT{//(fc,)FT[/(x)]} 


( 4 ) 
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This relation can be easily obtained from the well-known formula, FT{FT[g(x)]} OC g(-x). Eq ?? 
implies that the output will be proportional to the mirror image of the desired output function 

g{x)^ 

For performing Fourier transform we use graded index (GRIN) lenses. Since the optical prop¬ 
erties of dielectric GRIN lenses change gradually, the scattering of the wave could be significantly 
reduced which leads to higher efficiency. To realize the graphene-based type of such lenses the sur¬ 
face conductivity of graphene should be properly patterned in a way that the effective mode index 
of GP waves follows the quadratic refractive index distribution of their dielectric counterparts.^^ 

GRAPHENE-BASED META-TRANSMITT-ARRAY 

In this work, we implement the appropriate transfer function, H(kx) by means of a new type of 
meta-transmit-arrays on graphene. In order to manipulate the transmitted wave efficiently, not 
only the transmission amplitude should be completely controlled in the range of 0 to iJ ^^F^ Ifint 
also transmission phase should cover the whole In range independently.l^^^^ To this end, simi¬ 
lar to the approach of Monticone et a\.^ a meta-transmit-array comprised of symmetric stack of 
three metalines separated by a quarter-guided wavelength transmission line is utilized (Figure |^) 
in which each unit cell operates as a nanoscale spatial light modulator. To simplify the design 
procedure, the two outer stacks are chosen identical, but different from the inner one. The quarter- 
guided wavelength spacing between these stacks ensures that the transmission amplitude ripples 
are kept mi nimi z ed.!^ 

In order to fully control the transmission phase, in addition to amplitude, we need to locally 
manipulate propagating GP waves along and across the meta-transmit-arrayAs described pre¬ 
viously, this GP surface wave engineering is achieved via surface conductivity variation through 
an uneven ground plane beneath the graphene layer. 

Recently, an analytical results for the reflection and transmission coefficients of GP waves at 
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one dimensional surface conductivity discontinuity has been reported:^ 
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are the GP wavenumbers of left and right side regions of the discontinuity in the quasi-static ap¬ 
proximation, while Or^r represent their corresponding complex surface conductivities, and Eg is the 
average permittivity of the upper and lower mediums surrounding the graphene sheet. Similarly, 
for a GP wave incident on the discontinuity from the right side region, these coefficients, trr ans 
tRR, can be simply achieved by exchanging ki and kR in the eqs ??, ??. 

To relate the forward-backward fields on one side of the interface to those on the other side, 
the matching matrix is applied. By employing reciprocity theorem, the matching matrix of the 
interface is achieved as follows:!^ 


— 


^LR,m ^RL^m 
^LR^m ^RL,m j 


( 8 ) 


where m=0,.. .,5 is the interface of the building block (see Figure]^). Furthermore, the prop¬ 
agative forward-backward fields relation along a segment are defined by propagation matrix 


Pm = 




ikfjilm Q 
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where m=l,...,5 is the segment of the unit cell. Finally, the scattering parameters can be 
easily calculated using the whole building block transfer matrix obtained by multiplication of the 
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matching and propagation matrices 
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T=MQWPmMm ( 10 ) 

m=l 

The amplitude and phase of ^21 versus the chemical potentials of the inner and outer metalines, 
fic^in and lie,out ^ are plotted in Figure calculated by means of our analytical approach. It is 
obvious that by local tuning of the chemical potential of each unit cell any transmission phase and 
amplitude profile can be achieved. 

IMPLEMENTAION OF MATHEMATICAL OPERATORS AND 
REPRESENTATIVE RESULTS 

In this section, we implement first differentiation, second differentiation and integration operators, 
with the proposed structure. It is well-known that the derivative of a function is related to its 
first Fourier transform by: 

^^^ = F->{(aj“F(/W)} (11) 

Comparing eq ?? and eq ?? clarifies that to realize the derivation we have to perform a transfer 
function of {ikxY- So as described in the previous section we set the transfer function of the meta- 
transmit-array to H{x) oz {jxY. Since metalines are inherently passive media, the desired transfer 
function has to be normalized to the lateral limit to ensure that across the structure the maximum 
transmittance is unity, thus the appropriate transfer function is H{x )« (ix/iW/l))". Figures |^,c 
indicate the desired magnitude and phase profile of the first-order derivative transfer function /f(x). 
According to Figure by properly tailoring the values of the chemical potential along the lateral 
dimension for internal and external metalines, the transverse amplitude and phase distribution of 
the transfer function are implemented. To this end, the transverse distribution of chemical poten¬ 
tial for the designed first-order differentiator meta-transmit-array and its corresponding complex 
surface conductivity profile are depicted in Figure Now, a TM-polarized GP surface wave is 
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launched toward the designed meta-transmit-array and the calculated electric field distribution is 
shown in Figure As depicted in Figures |^,c the output profile of meta-transmit-array is in an 
excellent agreement with the desired transfer function. To be more precise, the standard deviation 
from the amplitude and phase of the desired transfer function is 0.04 and 6°, respectively. This is 
smaller than the standard deviation observed in Fig. S7b of which is about 0.19 and 15° for the 
amplitude and phase, respectively. As a result, a more accurate output response can be anticipated 
for the whole structures. 

Other operators can be designed in a similar manner. In order to perform the second-order spa¬ 
tial derivative, the desired transfer function is if (x) oc —{x/iW/2))'^. Although here the amplitude 
is a quadratic function of transverse dimension, the phase is constant. The results for the designed 
second-order differentiator are illustrated in Figure]^ 

To realize a second-order integrator, a challenge should be overcome. In this case, the desirable 
transfer function should be H{x) oc which leads to an amplitude profile with values tending 

to infinity in the vicinity of x = 0. To overcome this problem, we use the following approximate 
transfer function:^ 

{ 1 , if |x| < 

( 12 ) 

(/x)“^, if |x| > h 

where h is an arbitrary parameter that we set it to = VF/12.^ For all the points within |x| < h 
the amplitude is assumed to be unity, others follow the correct transfer function profile precisely. 
Figure shows numerical simulation of electric field distribution for this case. By comparing 
the obtained results from meta-transmit-array and analytical solution in Figures [^,c , an excellent 
agreement is observed. 

To demonstrate the functionality of our proposed structures, a GP surface wave in the form of 
a Sine function is sourced into our proposed GRIN/meta-transmit-array/GRIN configuration and 
the simulated electric field distributions are illustrated in figures|9lp^andpT]for the designed first- 
order derivative, second-order derivative and second order integrator, respectively. It is obvious 
from these figures that the achieved results are closely proportional to the desired results calculated 
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analytically. 


CONCLUSIONS AND OUTLOOK 

In summary, we have proposed and designed a new class of planar meta-transmit-array consisting 
of symmetric three stacked graphene-based metalines to perform wave-based analog computing. 
Using analytical results for the reflection and transmission coefficients of graphene plasmon waves 
at one dimensional surface conductivity discontinuity,we have demonstrated that full control 
over the transmission amplitude and phase can be achieved by appropriately tailoring of surface 
conductivity of each building block. Employing general concept of performing mathematical op¬ 
erations in spatial Fourier domain, we assign the meta-transmit-array a specific transfer function 
corresponding to the desired operation. The two designed operators in this work; i.e. differentiator 
and integrator, illustrate a high efficiency. The proposed graphene-based structure not only is ultra¬ 
compact, but also depicts more accurate responses than the bulky structure suggested in These 
features are due to exceptionally high confinement of surface plasmons propagating on a graphene 
sheet. However, this miniature size comes at a price, namely any future fabrication imperfections 
and tolerance will lead to distortion of the anticipated results, and therefore, degrade the efficiency 
of the proposed structure. The presented approach may broaden horizons to achieve more complex 
nanoscale signal processors. 

MATERIALS AND METHODS 

The full-wave simulations are carried out by the commercial electromagnetic solver, Ansoft’s 
HFSS. Since a high confined GP wave is stimulated here, in all our numerical simulations a free 
standing graphene layer is assumed.^ In the simulations, the graphene layer is treated as an surface 
impedance = 1 / cr^, where <7^ is the complex surface conductivity of graphene can be retrieved 
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by the Kubo’s formula under random phase approximation 
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where (O is the angular frequency, jlc is the chemical potential, T = l/2r is the relaxation time 
which represents loss mechanism, T is the temperature, e is the charge of an electron, h is the 
reduced Planck’s constant, and ks is Boltzmann’s constant. In our work, T = 300 K and T = 1 ps, 
which are relevant to the recently experimental achievements.^^ Moreover, A time dependence of 
the form is assumed and suppressed throughout this study. 

Due to very small mesh size, the simulation becomes computationally time consuming. As a 
result, we follow the approach in,l^i.e. block-by-block simulation. Initially, the first GRIN lens is 
excited with the desired input field profile. Then, the meta-transmit-array is excited by the moni¬ 
tored output field distribution of the first block. At last, the obtained output of the MTA is sourced 
into the second GRIN lens. To excite the GP surface wave with the desired field distribution, a 
current sheet is placed perpendicular to the graphene layer. By discretizing the current sheet to 
sufficient number of segments and assigning desired surface current amplitude and phase to them, 
each segment acts as an independent source. Consequently, any exact field distribution can be 
achieved easily. 

The schematic model used in the simulations is shown in Figure which contains two GRIN 
lenses and the meta-transmit-array embedded between them. The whole structure is restricted 
in transversal and longitudinal directions by W and 2Lg -\- D, respectively. Here, we set W = 
684 nm, Lg = 1028 nm, D = 100 nm, A = 18 nm, d = 5 nm. In all simulated examples, the 
wavelength is chosen A =6 jim. The proposed configuration can principally be implemented 
for other wavelengths using well-designed constituent parameters and dimensions, as long as the 
graphene layer supports the GP waves at that wavelength. A y-polarized GP surface wave with 
spatial variation f{x) = sinc(16;rx/IT) is used as the input function to evaluate the efficiency of 
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the whole GRIN Lens/MTA/GRIN Lens structure. |5] 
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Figure 1: Perspective view of wave propagation along a well-designed GRIN Lens/MTA/GRIN 
Lens which performs optical analog computing. 



Figure 2: Sketch of linear transversely invariant system to perform mathematical operations which 
consists of spatial Fourier transforming (FT), desired transfer function and inverse Fourier trans¬ 
forming (IFT) blocks. 
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(a) GRIN Lens MTA GRIN Lens 



Figure 3: (a) Schematic of two dimensional graphene-based computing system including GRIN 
Lens/MTA/GRIN Lens together with the input and output functions which variate across x direc¬ 
tion and propagate along z direction, (b) Basic building block of the metalines labeled by the 
propagation and matching matrices corresponding to the interfaces and segments, (c) Sketch of 
the meta-transmit-array made of three symmetric stacked metalines embedded between two GRIN 
Lenses. 



Figure 4: (a) Phase and (b) amplitude of transmission coefficient versus the internal and external 
metalines’ chemical potentials calculated by the proposed analytical approach. 
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Figure 5: The (a) transverse distribution of chemical potential for the designed first-order differen¬ 
tiator meta-transmit-array and (b) corresponding complex surface conductivity. 
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Figure 6: (a) Snapshot of electric field distribution {E^ for a TM-polarized GP surface wave in¬ 
cident on the designed first-order differentiator meta-transmit-array. Comparison of the transverse 
(b) amplitude and (c) phase distribution of transmitted wave and desired transfer function response 
right behind the structure. 
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Figure 7: (a) Snapshot of electric field distribution (E^) for a TM-polarized GP surface wave inci¬ 
dent on the designed second-order differentiator meta-transmit-array. Comparison of the transverse 
(b) amplitude and (c) phase distribution of transmitted wave and desired transfer function response 
right behind the structure. 
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(b) Transverse Amplitude Profile (c) Transverse Phase Profile 



Figure 8: (a) Snapshot of electric field distribution (E^) for a TM-polarized GP surface wave in¬ 
cident on the designed second-order integrator meta-transmit-array. Comparison of the transverse 
(b) amplitude and (c) phase distribution of transmitted wave and desired transfer function response 
right behind the structure. 



19 



























































1 




(D 


N 


a 

o 



N 

w 

HH 


0) 

N 


;-i 

O 

:z; 


Output real part 


0 


-1 





Figure 9: The (a) real and (b) imaginary parts of the output electric field compared with the ana¬ 
lytical results for the first-order differentiator, (c) Snapshot of the z-component of the electric field 
distribution along the GRIN Lens/MTA/GRIN Lens. The input function is f{x) = sinc(16;rx/VF). 
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Figure 10: The (a) real and (b) imaginary parts of the output electric field compared with the analyt¬ 
ical results for the second-order differentiator, (c) Snapshot of the z-component of the electric field 
distribution along the GRIN Lens/MTA/GRIN Lens. The input function is f{x) = sinc(16;rx/VF). 
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Figure 11: The (a) real and (b) imaginary parts of the output electric field compared with the ana¬ 
lytical results for the second-order integrator, (c) Snapshot of the z-component of the electric field 
distribution along the GRIN Lens/MTA/GRIN Lens. The input function is f{x) = sinc(16;rx/VF). 
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